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The stochastic response of nanoscale osciUators of arbitrary geometry immersed in a viscous fluid 
is studied. Using the fluctuation-dissipation theorem it is shown that deterministic calculations of 
the governing fluid and solid equations can be used in a straightforward manner to directly calculate 
the stochastic response that would be measured in experiment. We use this approach to investigate 
the fluid coupled motion of single and multiple cantilevers with experimentally motivated geometries. 
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Single molecule force spectroscopy using nanoscale 
cantilevers immersed in fluid is a tantalizing experimental 
possibility P, 01 • The precise manner in which nanome- 
chanical devices will be utilized for single-molecule force 
spectroscopy and sensing is currently under develop- 
ment l^ljl, however the detection system will rely upon 
the change in response of the cantilever due to the bind- 
ing of target biomolecules. It is therefore important to 
build a baseline understanding of the motion of a fluid 
loaded cantilever, or arrays of cantilevers, in the absence 
of active molecules. This is the focus of this Letter. 

The dynamics of the nanoscale structures considered 
are dominated by Brownian fluctuations although the 
mechanical structures are still large compared to the 
molecular size of the fluid molecules. The elastic re- 
sponse of a single cantilever immersed in fluid has been 
investigated previously in the context of atomic force mi- 
croscopy (see for example 0, 0)- This theoretical ap- 
proach, essentially two-dimensional in nature, modelled 
the cantilever as an inflnite cylinder oscillating in an 
unbounded fluid. This approach has been experimen- 
tally validated for long and slender micron-scale can- 
tilevers j^. However, the response of nanoscale can- 
tilevers (very strong fluidic damping) or short and wide 
cantilevers (where end effects would be important) is not 
well understood. Using the fluctuation-dissipation theo- 
rem we show that deterministic numerical simulations of 
the cantilever-fluid system can be used to obtain experi- 
mentally relevant stochastic quantities such as the noise 
spectrum of the fluctuations. 

The fluctuation-dissipation theorem allows for the cal- 
culation of the equilibrium fluctuations using standard 
deterministic numerical methods. This is possible be- 
cause the same molecular processes are responsible for 
the dissipation and the fluctuations. In the case under 
consideration here, predominantly these are the collisions 
of the fluid molecules with the cantilever, although dis- 
sipative processes within the elastic material of the can- 
tilever could also be included. 
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The fluctuation-dissipation theorem comes in many 
forms. For the fluid-damped motion of nanoscale can- 
tilevers it is sufficient to use a classical formulation. The 
most convenient form is the one originally discussed by 
Callen and Greene (see also [s^). Consider a dynam- 
ical variable A. For the classical system of interest here 
this will be a function of the microscopic phase space 
variables consisting of 3A^ coordinates and conjugate mo- 
menta of the cantilever, where N is the number of parti- 
cles in the cantilever. We investigate the situation where 
a force f{t) that couples to A is imposed. In this case 
the Hamiltonian of the system is 



H^Ho-fA 



(1) 



and we look at the linear response for very small /. Then 
it can be shown for the special case of a step function 
force given by. 



. I F for t < 
' ~ 1 for i > 



(2) 



that in the linear response regime, the change in the av- 
erage value of a second dynamical quantity B (again any 
function of the 3A'^ coordinates and momenta) from its 
equilibrium value in the absence of / is given by 



A{B{t))^PF{SAiO)SB{t)), 



(3) 



where /3 = {ksT)^^, ks is Boltzmann's constant, T is 
the absolute temperature, 6 A = A — {A)^, 6B — B — 
(_B)q and the subscript zero on the average () denotes the 
equilibrium average in the absence of the force /. Thus 
we can calculate a general equilibrium cross-correlation 
function in terms of the linear response as 



{5A{0)6B{t))^ = keT 



F 



(4) 
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There are no approximations involved in this result, ex- 
cept that of assuming classical mechanics and linear be- 
havior. If in addition the dynamical variables are suf- 
ficiently macroscopic that the mean {B{t)) can be cal- 
culated using deterministic, macroscopic equations, we 
have our desired result. 
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There are a number of advantages to a formulation of 
the fluctuation-dissipatfon theorem in time rather than 
frequency for our purposes. First, the fuU correlation 
function is given by a single (numerical) calculation, the 
response to removing a step force. The spectral proper- 
ties can be obtained by Fourier transform. This is partic- 
ularly advantageous for the low-quality factor (Q) situa- 
tion characteristic of fluid-damped nanoscale cantilevers, 
since the spectral response is very broad, and a large 
number of fixed-frequency simulations would be needed 
to characterize this response. A second advantage is that 
no expansion in modes of the oscillator is needed. Al- 
though such an expansion is not too hard for a high-Q 
oscillator where the dissipation has negligible influence 
on the mode shape, for small cantilevers in a fluid, the 
coupling to the fluid is large, and the motion of the fluid 
complex, so that a mode analysis would be quite difficult. 
Finally, the expression Eq. I^J allows us to calculate the 
correlation function and noise spectrum of precisely the 
quantity measured in experiment, firstly by tailoring the 
applied force to couple to one physical variable measured 
in the experiment (A), and then determining the effect 
on the second physical variable (B). This idea has been 
exploited in the very high-Q situation of the oscillators 
used in gravitational wave detectors although there 
it was convenient to formulate the result in the frequency 
domain. In our case for example, if the displacement of 
the cantilever is measured through the strain of a piezore- 
sistive layer near the pivot point 0, 3 of the cantilever, 
it is possible to tailor the force F to couple to this dis- 
tortion, and so determine the "strain-strain" correlation 
function of one or more cantilevers. 

Our scheme consists of the following steps in a deter- 
ministic simulation: (i) apply the appropriate force /, 
constant in time, small enough so that the response re- 
mains linear, and tailored to couple to the variable of 
interest A, and allow the system to come to steady state; 
(ii) turn off the force at a time we label t = 0; (iii) mea- 
sure some dynamical variable B{t) (which might be the 
same as A to yield an autocorrelation function, or dif- 
ferent) to yield the correlation function of the equilib- 
rium fluctuations via Eq. Q. Using the sophisticated 
numerical tools developed for such calculations, we can 
find accurate results for realistic experimental geome- 
tries that may be quite complex, for example the trian- 
gular cantilever design often used in commercial atomic 
force microscopy, or the reduced stiffness geometries cur- 
rently under investigation for use as detectors of single 
biomolecules as shown if Fig. ^ 

For simplicity, we illustrate our approach by finding 
the auto- and cross-correlation functions for the displace- 
ments Xi (i) of the tips of one or two nanoscale cantilevers 
with experimentally realistic geometries. To do this we 
calculate the deterministic response of the displacement 
of each tip, which we call Xi{t) i = 1,2, after switch- 
ing off at i = a small force applied to the tip of the 
first cantilever, Fi, given by Eq. For this case the 
equilibrium auto- and cross-correlation functions for the 
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FIG. 1: Schematic of the cantilever geometry (not drawn to 
scale): I = 3/^, w — lOOnm, h — 0.6/i, b = 33nm. The 
cantilever is silicon with a density ps — 2330Kg/m"^, Young's 
modulus Ea — 125GPa. This type of geometry is used to 
reduce the spring constant, the equivalent spring constant 
k = 8.7mN/m P. 

fluctuations xi and X2 are precisely 

(xi(t)xi(O)) = kBTXiit)/Fi, (5) 
(X2(i)xi(0)) = kBTX2it)/Fi. (6) 

The cosine Fourier transform of the auto and cross- 
correlation functions yield the noise spectrum Gii(i/) and 
Gi2{v), respectively, which are the experimentally rele- 
vant quantities. 

We have performed three-dimensional time-dependent 
numerical simulations of the deterministic fluid-solid cou- 
pled problem (algorithm described elsewhere 10, llj). 
The fluid motion is calculated using the incompressible 
Navier-Stokes equations and the dynamics of the solid 
structures are computed from the standard equations of 
elasticity. 

For long (/ > 100/i) and slender {I >> w) cantilevers, 
the cantilever response can be well approximated as an 
infinitely long oscillating cylinder 5, 6]. We first validate 
our numerical approach by investigating a cantilever in 
this regime. We emphasize that for the experimentally 
motivated nanoscale cantilevers of interest here an ap- 
proximate theory is not available, yet our numerical ap- 
proach remains valid providing a means to gain valuable 
insight. 

The micron-scale cantilever used for validation has the 
simple beam geometry as shown in Fig. |21 (b) (see case c2 
in 6]). For micron and nano-scale cantilevers immersed 
in fiuid, dissipation is dominated by the viscous motion of 
the fiuid driven by the cantilever vibrations. This can be 
described by a Reynolds number based on the frequency 
of oscillation ui as R — pwuP' /Arj where p is the fluid 
density and -q is the kinematic viscosity. For the can- 
tilevers of interest here the Reynolds numbers are typ- 
ically 0.01 ^ -R ^ 1 indicating that this is in the low 
Reynolds number regime. Small R corresponds to strong 
dissipation. 

The noise spectrum, Gii(i/), is calculated from the nu- 
merical results by taking the cosine Fourier transform 
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of the autocorrelations and is shown by the sohd line in 
Fig-EIa). The two broad peaks can be identified with the 
first two modes of the cantilever. The noise spectrum is 
also calculated using the long cylinder analytic approxi- 
mation and is shown by the dashed line. The analytical 
result for the fundamental mode of the noise spectrum is 
found in the following manner ^ (note that higher har- 
monics could be included if desired). In Fourier space 
the equation of motion for the cantilever displacement is 



k)x^Ff+FB, 



(7) 



where Ff is the force felt by the cantilever due to the 
fluid, 



Ff = meUj'^T{Lij)x, 



(8) 



Fb is the fluctuating (Brownian) force, TOe is the effec- 
tive mass 14i| of a fluid cylinder of radius w/2, m is the 
effective mass of the cantilever, and r(w) is the hydrody- 
namic function (for an infinitely long cylinder of diameter 
w oscillating in the x-direction) which contains both the 
fluid damping and fluid loading components and is given 

by El, 



r(w) = 1 



AiKi{-iViR) 



iRKoi-WiR) 



(9) 



where Ki, Kq are Bessel functions and i — \/— T. From 
the fluctuation-dissipation theorem the spectral density 
of the fluctuating force, Gpglv), can be related to the 
dissipation due to the fluid and is given by 



Gfb{v) = ■AkBTmeTQUjTi{uj), 



(10) 



where oj = 2ttv. Solving for the spectral density of 
the displacement fluctuations, Gnc{v), from Eqs. 0, © 
and yields. 



AkBT 1 



k 



UJQ 



(11) 



(1 - w2 (1 + ToTriRoCdW + {u^ToT,{Ro^)) 



where lj = lu/ujo is the frequency relative to the vacuum 
resonance frequency too = \/k/rn, Rq is the Reynolds 
number based on wo, To is the ratio of the mass of fluid 
contained in a cylindrical volume of radius w/2 to the 
mass of the cantilever and Fr and F, are the real and 
imaginary parts of F, respectively Sader's analysis 
does not take into account the frequency dependence of 
the noise force and assumes that the numerator is con- 
stant. The frequency dependence is not large for R < 1; 
however the correction has been included in our analysis. 
Equation yields the analytical curve for Gii(z^) in 
Fig. [2Ia) and the agreement is excellent with our numer- 
ical results. 
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FIG. 2: Panel (a) the noise spectrum from simulation (solid 
line) and from theory (dashed line) . For the theory calculation 
only the fundamental mode has been considered. Panel (b) is 
a schematic of the micron-scale cantilever used for validation: 
length I — 197/i, width w = 29/i and height h — 2/i. The 
applied step force is Fi — 26 nN. 
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FIG. 3: Schematic showing various cantilever configurations. 
In all configurations the step force Fi is released at t = 
resulting in the cantilever motion referred to by Xi{t). 
The motion of the neighboring cantilever is, X2{t), and is 
driven through the response of the fiuid. Panel (a) two can- 
tilevers with ends facing, panel (b) side-by-side cantilevers 
and panel (c) cantilevers separated along the direction of the 
oscillations 



A recent study of the fluid coupled motion of two adja- 
cent Ifj, beads illustrates the importance of understanding 
the cross correlations in the fluctuations for use in sin- 
gle molecule force measurements 13]. The fluid distur- 
bance caused by an oscillating cantilever is long range, 
producing motion of the other cantilevers through the 
viscous drag. As a result, the stochastic motion of mul- 
tiple cantilevers will be correlated. However, the numer- 
ical approach developed here remains valid for multiple 
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fluid-loaded cantilevers of arbitrary geometry, and our 
approach can be used to quantify the response of multi- 
ple cantilevers with the precise complex geometries used 
in experiment (as we show below) as well as to help de- 
velop a better analytical understanding of idealized ge- 
ometries. Various possible cantilever configurations are 
shown in Fig.|3Ia)-(c); here we will present results for the 
end to end case and defer results for the other geometries 
to a later paper. 




tilever fluctuations, given by the cosine Fourier transform 
of the cross-correlation function, are shown in Fig. EJa) 
and (b). Notice that tuning the separation could be used 
to reduce the correlated noise in some chosen frequency 
band. 
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FIG. 4: Absolute predictions of the auto- and cross- 
correlation functions of the equilibrium fluctuations of the 
cantilevers shown in Figs. Q and El Panel (a) a schematic of 
the step force boundary condition that is applied to the tip 
of the first cantilever, Fi — 75pN. Panel (b) shows the auto- 
correlation and panel (c) the cross- correlation of the fluctua- 
tions (5 separations are shown s — h, 2h, 3h, ih, 5h where only 
s = h,5h are labelled and the remaining curves lie between 
these values in sequential order). 

We use our approach to calculate the behavior of the 
experimentally motivated cantilever shown in Fig.^ Full 
three-dimensional simulations were performed for both 
one cantilever and two cantilevers facing end to end in 
fluid as shown in Fig. OJ^a). Through the fluctuation- 
dissipation theorem the simulations yield results for the 
cantilever autocorrelation function and the two cantilever 
cross-correlation function shown in Fig. ^b) and (c), re- 
spectively. The response of the driven cantilever (on the 
left) was not signiflcantly affected, for any of the separa- 
tions investigated, by the presence of the adjacent pas- 
sive cantilever on the right. The value of (a;i(0)a;i(0)) is 
0.471nm^ indicating that the deflection of the cantilever 
due to Brownian motion in an experiment would be 
0.686nm or about 2.3% of the thickness of the cantilever. 
The cross-correlation of the Brownian fluctuations of two 
cantilevers is small compared with the individual fluc- 
tuations. The largest magnitude of the of the cross- 
correlation is — 0.012nm^ for s = h and — 0.0029nm^ for 
s = bh. The noise spectra for both the one and two can- 



FIG. 5: Panel (a): The noise spectrum, Gi\{u). Panel (b): 
The noise spectrum, Guii'), as a function of cantilever sepa- 
ration, s, for two adjacent experimentally realistic cantilevers 
(see Fig. 0. Five separations are shown s = h, 2h, 3h, Ah, 5h 
where only s — h,5h are labelled and the remaining curves 
lie between these values in sequential order. 

The variation in the cross-correlation behavior with 
cantilever separation as shown in Figure ^c) can be un- 
derstood as an inertial effect resulting from the nonzero 
Reynolds number of the fluid flow, as follows. The flow 
around the cantilever can be separated into a potential 
component, which is long range and propagates instan- 
taneously in the incompressible fluid approximation, and 
a vorticity containing component that propagates diffu- 
sively with diffusion constant given by the kinematic vis- 
cosity u. For step forcing, it takes a time Ty = s^/j^ for 
the vorticity to reach distance s. For small cantilever 
separations the viscous component dominates, for nearly 
all times, and results in the anticorrelated response of the 
adjacent cantilever in agreement with (ISjj. However, as 
s increases the amount of time where the adjacent can- 
tilever is only subject to the potential flow fleld increases 
resulting in the initial correlated behavior. 

In summary, a numerical approach to calculate the 
stochastic response of single and multiple nanocan- 
tilevers, of arbitrary geometry, immersed in a viscous 
fluid using deterministic calculations has been developed, 
validated, and applied to an experimentally relevant can- 
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tilever geometry. The methods described here are apph- 
cable to atomic force microscopy in general and also to 
other nano-structure fluid interaction problems which are 
of growing importance as NEMS technology advances. 
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